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Abstract. We study energy measures on SG based on harmonic functions. 
We characterize the positive energy measures through studying the bounds of 
Radon-Nikodym derivatives with respect to the Kusuoka measure. We prove 
a limited continuity of the derivative on the graph V* and express the average 
value of the derivative on a whole cell as a weighted average of the values on 
the boundary vertices. We also prove some characterizations and properties 
of the weights. 



1. INTRODUCTION 



The Sierpinski gasket, or simply SG, is the unique nonempty compact set satis- 
fying 

2 

SG = (J F.SG 

i=0 

where F t = ^(x + qi) and {<Zi}?=o are vertices of an equilateral triangle in counter- 
clockwise direction. {%}f =1 are called the boundary points of SG. If w = (wi, • • • , w n ) 
is a finite word, we define the mapping F w — F Wl o ■ ■ ■ o F Wn where \w\ is defined to 
be the length of the word. Define V a — {qt}, V n = \JFi (V n -i) and V* = U^Li ^n- 
For the boundary of a cell C — F W SG , denoted by dC, we mean the three vertices 
{F w (q l )}f =0 of the cell. 



For function u, v defined on SG, we define the energy on level m by 



(1.1) £ m {u,v)=[-) ^2 ^2(uoF w (q i )-uoF w (q j ))(voF w (q i )-voF w (q j )). 

\w\=m i<j 



and energy E(u, v) — lim £. m (u,v) provided that E(u,u) and E(v,v) are finite. 
For each cell F W SG, we define the energy of u, v on this cell by 

^F^sciu^) = l-j £.(uoF w ,voF w ). 
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Observe that the energy is non- negative if u = v, but in general it is not. In the 
case u = v, we write £(u) instead of £(u,u). We define a function h to be har- 



monic if it minimizes the energy from level m to level m + 1 , as defined in (1.1). A 
simple computation shows that this minimization occurs when £ m (/i) = £ m+ i(/i) 
at which to each cell F W SG, h(F w0 ( qi )) = lh{F w (q )) + ^h(F w ( qi )) + \h(F w {q 2 )). 
Equivalently, a function h on SG is harmonic if the sequence, taking u = h = v 
in (1.1), is a constant sequence, i.e. constantly ^2 i< j(h(q i ) — h(qj)) 2 . Given three 



values on the boundary, in order to have a constant sequence in ( |1.1[ ), all the other 
values are completely determined. Consequently, all the harmonic functions form a 
three-dimensional space and hence any values on the boundary of a cell can extend 
to a unique harmonic function on SG. A harmonic function h is said to be sym- 
metric(resp. skew-symmetric) if h(qi) — h(qi + i)(resp. h(qi) — —h(qi + i)) for some 
i. Harmonic functions characterized by hi(qj) = Sij are three symmetric functions 
and the notation hi will be used throughout this paper. On the space 5f of har- 
monic functions modulo constant, (JC, £) is a two-dimensional inner product space. 
A harmonic function h is said to be symmetric (resp. skew-symmetric) in 'K if there 
exists a constant c and a symmetric (resp. skew-symmetric) harmonic function / 
such that h = c + /. Explicitly, for a symmetric harmonic function h, we mean 
h = ahi + c for some a, c £ K and some and for a skew-symmetric harmonic 
function h, we mean ahi — bhj + c for some a, b being positive, c € R and some 
hi, hj. To each harmonic h, we write h 1 - the harmonic function having the same 
energy of h and orthogonal to h unless we specify h "orthonormal to h" in which 
case we mean having energy 1 and orthogonal to h. 

The standard measure /i is a self-similar measure characterized by [i(F w FiSG) = 
^fi(F w SG) for every cell C. The standard Laplacian Am of u is a continuous 
function satisfying 

(1.2) - l{u,v) = J Auvdn 

for every v with £(i>) < oo and vanishing at the boundary. Harmonic functions on 
SG are precisely those whose Laplacian is zero. However, for a function u whenever 
Au is defined as a function, A (it 2 ) is not. We can view Am as a measure m = Audfi 
and Au exists as a function if and only if this measure is absolutely continuous with 
respect to the standard measure. The carre du champs formula given by, for all 
/, u, v of finite energy, J SG fdv U:V = ^E(fu, v) + ^£(u, fv) — ij£(/, uv) shows that 

£(u 2 ,u) = —2 J (uAu)vd/i — 2 J vdv u 

for all u, v with Au as a function. Thus if Au exists as a function and if we first 
view A(u 2 ) as a measure, 

A(u 2 ) = 2uAudfi + 2dv u 

and A(u 2 ) exists as a function if this measure is absolutely continuous with respect 
to the standard measure /i but this is not the case. 

Suppose that u, v are functions of finite energy. We can define a signed measure 
v u ,v by, on each cell, 

Vu,v(C) = £. c (u,v). 
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As usual, if it happens that u — v, we write v u instead. In this paper, we will 
restrict the attention to both u, v being harmonic. The energy measures of harmonic 
functions form a three-dimensional space. The Kusuoka measure v is defined to be 
+ v\ + v%. By [5J, and by some normalization factor, we know that + v h ± = cv 
where c = |£(/i). Now, we define the energy Laplacian A„u of u by, for every finite 
energy v vanishing on the boundary, 

(1.3) - £(u,v) = J (A u u)vdv. 

By [5J, every energy measure of harmonic functions is absolute continuous with 
respect to the Kusuoka measure v and it makes sense to consider the Radon- 
Nikodym derivatives. For the energy measures of the symmetric harmonic func- 
tions, we will denote by vi instead of i>h v \v{\ forms a basis of the space of 
all energy measures of harmonic functions. Since ho + h± + h% = 1, we have 
Vho,hi + Vh u hj + Vh 2 ,hj = v hhi = for j = 0, 1, 2. Since v^ g = u g j for all /, g, this 
gives us three equations that we solve to obtain 

We will keep using the notation v^ g for any energy measure of harmonic functions 
without mentioning / and g. 

In the development of analysis on fractals, three basic concepts are energy, mea- 
sure and Laplacian. For SG, which may be regarded as the "poster child" for the 
theory, Kigami has developed an elegant construction of a self-similar Laplacian 
based on a self-similar energy £ and self-similar measure /i via an analog of integra- 
tion by parts. This theory is briefly sketched below, and is described in complete 
detail in the book [8]. This standard Laplacian A is not the only Laplacian on 
SG that has been widely studied. A competing energy Laplacian A„, based on 
the same energy but using the Kusuoka measure v in place of fi, has received a 
lot of attention in recent years([H |5j El [101 E]). It has two striking advantages 
over A, namely the existence of a Leibniz-type product formula and the fact that 
the associated heat kernel satisfies Gaussian type estimates. The main drawback 
so far has been that the Kusuoka measure is not algorithmically transparent. The 
purpose of this paper is to eliminate this drawback. We will provide algorithms 
that describe the nature of this measure and more generally the family of energy 
measures based on harmonic functions h, and the Ra don- Nikodym derivatives 



Some of the results of this paper, especially Theorem 5.5 may also be derived 
from [B]. 

By pp, we have that there exist matrices Ei such that for every cell C, 

v f , g (F t2 C) J \ v Lg {F 2 C) 

The EiS are all diagonalizable and so easy to calculate E™. Explicitly, after diag- 
onalization, we have 
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Eq = I I 1 



3 \ / A 



3 
5 

ii i / V o o 



(1.5) E 1 = i 3 



£ 2 =l 7 

11 / V 
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By symmetry of SG, it is easy to calculate that the equality of limits of quotients 
which are to obtain the value of derivatives 

(111) E^E 1 \ y ) (111) ETEo 

lim 



where 




(111) Eft 1 Ei \ \ (111) E™E 



a \ 








; - 





It immediately shows that it does not matter where the sequence of cell is taken. 
For each infinite word w, denote w m be word of the first m letters of w, then 
lim E w = 0. It follows that every energy measure is continuous. If we renor- 

m— >oo m 

malize E w by (I)' 1 "', which is inverse of the largest eigenvalue of E^, then we get a 
non-zero limit for each After computation, the limits are: 



in, [ I ) E, 







(1.6) lim (|) E? 



>3C 



2 



42 


-3 


-3 


14 


-1 


-1 


14 


-1 


-1 


-1 


14 


-1 


-3 


42 


-3 


-1 


14 


-1 


-1 


-1 


14 


-1 


-1 


14 


-3 


-3 


42 



Notice that all of the above are rank-one matrices. 

It is natural to ask whether the energy Laplacian behaves similar to the standard 
Laplacian; in Section 2, we express the "self-similarity" of the energy Laplacian. 
In Section 3, we will study the decay rates of energy measure from cells to its 
subcells. Through studying this, we give sharp bounds for the Radon-Nikodym 
derivatives for each v^. Being in general a signed measure, it is natural to ask when 
it is a (positive) measure; in Section 4, we will characterize all the positive energy 
measure. It is well-known that the derivative is not continuous; nevertheless, in 
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Section 5, we find out some restriction of the derivative is continuous. In Section 6, 
we will express, on each cell, the average value of derivatives by a weighted average 
of values on vertices of the cell. How the energy distributes is mysterious; in Section 
7, we will provide some graphs and hypothesis about these. 



2. Self-similarity of energy Laplacian 



In this section we discuss the self-similarity identities for energy measures and 
the energy Laplacian. We will see that the Radon-Nikodym derivatives Ri = ^ 
play a crucial role in these identities. Individually, none of the Vi is self-similar; 
nevertheless, together they form a self-similar family, in the sense of Mauldin- 
Williams[7]. We begin with considering the symmetric functions ho- We have the 
following equalities, 



(2.1) 



2 2 3 
h oF = h + -hi + -h 2 — 1 + -h a ; 

2 1 
h Q oF 1 = -ho + -h 2 ; 

2 1 
h o F 2 = -h + -h x . 



Next we establish the relation of i/j's from cells to subcells. 



Theorem 2.1. Let 



x / 9 \ / 2 2 -1 \ / 2 -1 2 \ 

M = - 2 2 -1 ,Mi = — 9 ,M 2 = - -1 2 2 

10 \2-12/ 10 \ -122/ 10 \ 9 / 



then 



vi{FiC) \=M i \ Vl {C) 
v 2 {FiC) J V "2(C) J 



for every cell C. In other words, 



v 2 



= Mo 



Vo 
V\ 

v 2 



o Fq 1 + Mi 



Vo 
V\ 

v 2 



o Ff 1 + M 2 



Vo 
V\ 

v 2 



o F 
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Proof. Let / be a contintuous function on SG, by (1.4 1, (2.1), 
fdv Q 



F SG 



SG 



f o F Q du ho+ 2 hl+ 2 h2 



I ° F du 



SG 



25 



/ o F dv x + — 

SG 20 



/ o F Q dv 2 



SG 



2 

+ 5 



SG 

4 

/ o F Q dv + — 

SG 1)0 

f o F dv - 

SG 



f o F Q dv hoM + - fo F dv hoM + 



SG 

4 

/ o F dvi + — 

SG 20 



25 



fo F du huhl 



SG 



f o F Q dv x 



SG 



SG 



f o F Q dv 2 

SG 

f o F dv 2 ) 



SG 



+ ~ ( - / / o F Q dv + fo F dvx - f o F du 2 ) 



SG 



SG 



25 



SG 



f o F dvi 



SG 



f o F dv 2 



SG 



f o F a dv . 



I SG 

On the other hand, 



fdv 



F\SG 



5 
3 

5 / 1 



SG 



f o Fidvi ho+ i h2 



3\25j SG J 25 



SG 



f oF x dv 2 + — 



2 

15 



/ o Fidi'o + 



S'G 



15 



/ o Fidz/i - 



SG 



15 



/ o F 1 di/ ho ,h. 

SG 

f o Fidi/ 2 



SG 



using (1.4) again. Similarly, 
fdvQ 



F 2 SG 



5 
3 
2 

15 



/ o F 2 dvi ho+ i h2 



SG 



f o F 2 dv Q - — 

SG 15 



/ O F 2 ^l + T7 
SG it( 



/ o F 2 dv 2 . 



SG 



Since this is true for arbitrary continuous function /, we have the following identity 

9 



15 
2 

15 
2 



u o F 
i/ o Ff 



2 , 1 , 

15 1 15 1 



-i 1 -i 2 

VnOF n l/i of, H V? O F 

15 2 15 2 15 2 
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By symmetry, we get the similar expressions for v\ and vi , and we have the relation 

vo \ ( "a \ ( \ I v Q 

vx =M vi of^ + Mi 1/1 \oF^ + M 2 \ vi I oFj" 1 . 



f 2 / \ v % \ v 2 \ V% 



□ 



Corollary 2.2. TTie Kusuoka measure satisfies the variable weight self-similar 
identity 



2 ' ' 1 12 



i=0 



Proof. If we consider j/ = ( 1 1 1)| ^i | • where i?j = the Radon- 



;/ 2 



Nikodym Derivative of Vj w.r.t. we get 
(2.3) 

12 -i 1 i 12 , 1 i 12 , 1 

v = — vq o F n 1 H 1/ o F n H j/i o F-i H f o F H z/ 2 ° F> H f ° F> 

15 15 15 1 15 1 15 2 15 2 

□ 

We now give the "self-similarity" of the energy Laplacian A y . 



Theorem 2.3. Let Q 3 , = ± + ^Rj. Then the energy Laplacian A„ satisfies 

A v {uoFj) = Qj(A v u) oFj 

Proof. For any integrable /, we can write 



j=o- 



which is equivalent to 



fdv= j (^ + ^ R \foF j dv. 

F 3 SG JSG 10 / 

By this relation and the self-similarity of v, we can have the self-similarity of A u u, 
that is to relate A v (uo Fj) and (A v ii) o Fj. By definition of energy Laplacian and 

by g, 

—E(u,v)= / (A u u)vdv 

JSG 



J=Q JSG it( 15 
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On the other hand, the self-similarity of £ also gives 

5 2 

-E(u,v) = --^ZiuoF^voF,). 

Together we have 

2 2 

/ J2(t^ + T5 R 3)( A - u ) oF j voFj du=^Y^f A^uoFj) voFj dv. 

JsG j=0 j=0 JsG 

Since v is arbitary, and so v o Fj is arbitary, we must have the self-similarity of A v 

3 

-Qj{^ v u)oFj = A^uoFj). 
o 

□ 

Corollary 2.4. Let w = (wi, ...,w m ) be a finite word of length m and F w — 
F Wl o .... o F Wm . Define 

Q w = Q Wm ■ (Qw m -! ° F Wm ) ■ (Q Wm _ 2 o F Wm _ 1 o F Wm ) ■ ■ ■ (Q Wl o F W2 o • • • o F Wm ) 

Then 

A v (u o F w ) = Q w (A v u) o F w . 

Proof. Use the theorem and iterate. □ 

To compute A v u on the cell F W SG, zoom in by looking at u o F w as a function 
on SG, compute its Laplacian A v (u o F w ) and multiply by 

(A„it) oF w = ^-A v {u o F w ), 

and then zoom out via F~ x 

^uu\ FwS G = f^-A„(« o F w )\ o F- 1 = 1 (A„(u o F w )) o F-\ 
\Qw J Qw°r w 



3. Bounds of derivatives 

By the definition of energy measures, we see that every energy measure is con- 
tinuous; by the self similarity of SG, we see that it makes sense to consider the 
measure from cells to subcells. It is natural to study how the measure varies for a 
sequence of cells converging to a point. 

Lemma 3.1. For a harmonic function h, if we let x = v h {F w F SG), y = v h (F w FiSG), 
z = VhiFwFzSG), then 14x — y — z > with equality if and only if h is skew sym- 
metric (with respect to F w q ) on the cell F W SG. 
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Proof. WLOG, let h(F w q ) = 1, h{F w qi) = a, h(F w q 2 ) — 0. Then, writing m = \w\ 
and a direct computation on the energies on each cell, 



y = 



5\ m+1 ( (> 
3, 

5\ m+1 / 6 



) (l)(° 2+ " + 1 > ; 

j (3a 2 - 3a + 1) ; 



3/ V 25 

m+l 



so lAx-y-z = 10(f)" 1 (4) (a 2 -4a + 4) anda 2 -4a + 4 = (a-2) 2 > 0. Equality 
holds iff a = 2, in which case h is skew symmetric. □ 

Lemma 3.2. Let w be a word. For 









■ 




i 









we /lave 

14x - y - z > 0. 

Proof. Since v — J^fj, by Lemma 3.1, we see that 14vi(F w F(f) — Vi(F w F\) — 
Vi(F w F2) > for all i. To complete the proof, we show that \Avi(F w F^) — 
Vi(F w F\) — Vi(F w F2) > for some i. Suppose, on the contrary, that each hi 
is skew symmetric with respect to F w qo on the cell F W SG. Then there exist Cj,aj 
such that (hi o F w (q ), hi o F w (q 1 ),h i o F™^)) = (cj,Cj + a*, Cj - a») = Cj(l, 1, 1) + 
aj(0, 1,— 1), but this shows hi lie in a two-dimensional subspace, contradicting hi 
being linearly independent. □ 

Lemma 3.3. For the sequence {F w F^ n SG} m , which converges to the point F w (qi), 
we have that v(F w FV rl SG) = 6 ((f)" 1 ). 

Proof. We prove this statement for i = 0; the other cases are similar. 
We know that 

v(F w F™SG) = v (F w F™SG) + Vl {F w F™SG) + v 2 (F w F™SG) 

6/5 

(111) E™E U " 



= ( 1 1 1 ) E™E„ 

Now we let 

( 2 
E w I 2 

V 2 

and see from the diagonalization of E™ that 

, { f^sg) .(-»,+ •, + •,) (1)"' + i (i4* - , - „ (I)" . e (I 




2/5 \ 


/ 2/5 


6/5 - 


- 2/5 


2/5 / 


V 6/5 
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since 14x — y — z > by the Lemma 3.2. □ 



We have proven the decay rates of the measures of a sequence of cells converging 
to a point which is not equivalent to a skew-symmetric cell. Surprisingly, in the 
case which the cell is equivalent to skew-symmetric at that point, we have the first 
condition for the Radon-Nikodym Derivatives being 0. 

Lemma 3.4. If h is symmetric on the cell F W SG with respect to the point F w qo, 
v h {F 2 F™SG) = e((^) m ) and hence ^(x) = if x = F w Ft(q 2 ). Conversely, 
if fa^{F w F\(q2)) = 0, h is symmetric on the cell F W SG with respect to the point 
FwQo- 



Proof. We know that Vh(F w SG) = (§) v hoF w , so we may consider the case when 
w is the empty word. We see from Lemma 3.1 that the term (|) m in Vh(FiF™) 
vanishes if and only if h is skew-symmetric on FiSG with respect to the point 
Fi_F , ™(g 2 ) = -Fife), which is true if and only if h is symmetric on SG with respect 
to the point qo, in which case Vh {FiF^SG) =Q((j E ) m ) and the denominator of 
^f(x) = lim ^(fyF™ sg) dominates. It follows that h is symmetric on the SG is 
sufficient and necessary to have ^(-^1(92)) =0. □ 



Theorem 3.5. let h be a harmonic function on SG with Vh = ovq + bv\ + cv 2 and 
h be a harmonic function orthonormal to h under the energy inner product. Also 
let C be a cell in SG. Then 
a) 

inf ^(*) = 0; 
xec dv 



sup — = -(a + b + c). 

xec av 3 



and if the maximum of ^f is attained, then the minimum of —^j- is attained at 
the same point. 



Proof. We first prove a). We know from j3] that for a harmonic function h with 
h(v ) = 0,h(vi) = l,/i(u 2 ) = a and a > 2 (WLOG, as argued in [3]), there is a 
continuous, monotone decreasing function M which maps a to the unique point x 
on the edge connecting Vq,Vi where h attains its maximum (relative to that edge). 
We view this edge as an interval, identifying v with and V\ with 1. 
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We also know that for the symmetric function h 2 , the value of the derivative at 
Fi(vq) is zero. We will now show that if a; is a dyadic point along the edge, the 
value of the derivative of h is zero. 

If a; is a dyadic point, then x £ V m for some m and so z ~ m x ~ m y such that 
y = F w (v ), z = F w (vi), and x = F w Fi(v ). 








When we restrict to F W SG, we have a harmonic function with maximum value (on 
the left edge) at F w (F 1 (q )). As discussed in [3], h 2 is the unique harmonic function 
such that the maximum is attained halfway along the edge, so the derivative at x 
is zero . 

Now suppose that x is not a dyadic point. 

Since dyadic points are dense in the interval [0,1], we pick a increasing sequence 
x n of dyadic points converging to x. Because M is decreasing and continuous, there 
exists a decreasing sequence of a n converging to a satisfying M{a n ) = x n . 

Let e > be given. A simple computation shows that the energy measure i>h = 
avQ + (1 — a)v\ + (a 2 — a)v2- From the previous discussion, we can pick an a n (x n , 
resp.), which will be denoted as a e {x e , resp.), satisfying \a— a £ | + |(l — a) — (1 — a e )| + 
I (a 2 — a) — (a 2 . — a £ )\ < e. We also denote the harmonic function with boundary 
values 0, l,a £ as h e . The previous discussion shows that d ^ L - — 0. Now, 
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dv 

■ d(Vh : 1/hJ (x £ )+ d, ^( Xs ) 

dv dv 



dv 

= (o - a e )^{x e ) + ((1 - a) - (1 - a e ))^{x e ) + ((a 2 - a) - (a 2 - a*))^(* e ) 
<|a - a e | + |(1 - a) - (1 - a e )| + |(a 2 - a) - (a 2 - a e )| 



where the penultimate inequality comes from the fact that < 1 since z/,(C) < 
i>((7) for any cell C . 
And consequently, 

inf p- = 0. 

On a general cell C, /j|c is harmonic, so the same argument applies. 

We then show b). By [8J, we can write |f = tv-h + v h ±, where r 2 = ^/gg) • Since 
the total measure of Vh on SG is 2 (a + 6 + c), we have 

3 c% , ^dv h ± 



2 (a + b + c) dv dv 
It follows that 



Whence, 



3 dv h . 

— : , r SUp — — = 1 — 6 mi : = 1. 

2(a + o + c) X £c dv xec dv 



dvh 2(a + 6 + c) 

sup — = . 

xeSG dv 6 



Further, if the derivative attains its suprcmum, then the maximum occurs where 
occurs in *V" , that is, where the local extremum occurs in the smallest edge of 
the harmonic function h . □ 



As we can see, on every cell of SG, the supremum and infimum are distance 
2(a+b+c) £ rom gg^h ther, so the derivative is far from continuous. However, we will 
later recover limited continuity when we restrict to the edges of cells. 



4. Characterzation of Positive Energy Measures 

Within the 3-dimcnsional space of signed energy measures which measures are 
positive? The following result gives the answer. 

Theorem 4.1. Let P be the set of {a, b, c) £ M 3 such that avg+bvi+cv2 is a positive 
energy measure. Then P is the solid cone S = {(x, y, z) g M 3 : xy + yz + xz > 0}, 
and the boundary of P is the set of (a, b, c) such that av + bv\ + cv 2 — v^ for some 
harmonic function h. 
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Proof. We first show that S C P. We know from pQ that the set of (x, y, z) such 
that x — Uh(F SG), y — v h {FiSG), z — v h (F 2 SG) form the cone 

I i | \ 2 2,2,2 

— (x + y + z) =x +y +z . 
A simple computation shows that 

transforms (x,y,z) to the corresponding coefficients (a, b, c). Hence, the set of 
(a, b, c) corresponds to some is the cone 

ab + be + ca = 0. 

Now we show that the coefficients that come from harmonic functions are on the 
boundary of P. Consider the energy measure Vh — ev for e > 0. Suppose, for 
contradiction, that this measure is positive. Then for any cell C, {y^ — ev)(C) > 0, 
so v ^c) > e > for all C, but we know the innmum should be zero, so this is a 
contradiction. Therefore (a — e, b — s, c — e) ^ P and (a, 6, c) is on the boundary. □ 

Immediately, since the solid cone is convex, every positive energy measure is pre- 
cisely a convex combination of two energy measures of a single harmonic function. 
But we can prove a little bit more. 

Corollary 4.2. Every positive energy measure of harmonic functions is a convex 
combination of positive energy measures of h and h 1 - for some harmonic function 
h, and vice versa. 

Proof. For every harmonic function h, Vh + Vh 1 - — cv f° r some c > 0. This shows 
Vh, Vf l ± , v ne on the same two-dimensional subspace. 

Suppose a two-dimensional subspace W containing v is given. W intersects the 
boundary of the cone S stated in the theorem on two and only two lines which 
contain energy measures of a single harmonic function. The above paragraph proves 
that one of the two lines contains Vh for some h while the other contains v h ± . 

Now, given any non-zero positive energy measure a (the zero case is trivial), con- 
sider the two-dimensional subspace V containing a and v. Let P be the plane 
passing through a normal to (1, 1, l)(in the sense of embedding the energy measure 
space into IR 3 ). And there are distinct functions h, f lying in the intersection be- 
tween P, V and S. Assume fh = Y] a i v i an d Vj — X) biVi- Since both /, g lie on P, 
E(h) = 2J2 a i — %J2bi — £(/)• And the above discussion shows / = h 1 - modulo 
constants. The converse is obviously true. □ 

5. Limited Continuity 

We have seen from the previous section that the derivative of an energy measure 
is not continuous; however, if we restrict the derivative to the set of vertices V*, it 
is continuous on the edges of every triangle. This can be seen from Figure |5.1| 
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1 



Figure 5.1. Radon-Nikodym Derivative of Uq with red dots indi- 
cating minimum values and green dot indicating maximum value. 



We know from ( 1.6 1 that 



40 I Hi), | | ] E'i 



( -1 -1 14 )E W A 



-1 


-1 14 \ 


-1 


-1 14 =: A 


-3 


™3 42 / 


~~ A m 


We also know 


c w ( 


-1 -1 14 ) 



( -1 -1 14 )E. 



and in fact 



We are going to prove, for any finite word w consisting only 1 and 2, that 
is bounded away from zero; to do this, we need the following lemmas. 




'5^1™! 



Lemma 5.1. Given any finite word w consists of letters 1 and 2 only, and letting 











= E W | 


■ 
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then 

' X \ I X 

( -1 -1 14 )E 2 [ y J > ( -1 -1 U)E 1 [ y 
Proof. We prove this by induction. It is trivial for w being the empty word. Let 

Suppose that the conclusion holds for \w\ — to, which is equivalent to saying x — 
Ay + llz > 0. Also assume for induction that z > x. Consider first the word wl. 
We have 

(-1 -1 U)(E 2 -E 1 )E 1 \ y J = | ( 1 -4 ll)E 1 ^y 

= ^(-9x-4y + 2lz) 

= ^((x - Ay + llz) + 10(z - x)) 
> 0. 

x 

Furthermore, the third entry of E\ | y | minus the first entry is 

15(0 -x) > 0. 

Now we turn to the second case concerning the word w2. 
( -1 -1 14 ) (E 2 - E X )E 2 I y 1=^(1 -4 11 ) E 2 I y 



z 



5 

l(-3y + 19s) 

^ ((x - 4y + llz) + (z - .t) + y + 7z) 



> 0. 

x\ 

Also, the third entry of i?2 ( V minus the first entry is 

^(-4x + y + llz) = — (y + 4(z - x) + 7z) > 0. 

The positivity of y and z comes from the fact that 

/ 1 \ / u 2 {F w0 SG) 
E w \ 1 = - u 2 (F wl SG) 
V 3 J 5 V V2(F w2 SG) 

This completes the induction. □ 
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Lemma 5.2. Under the same hypothesis and notations as in lemma 1, we have 



( -1 -1 14 )E 1 E 2 E W 1 > ( -1 -1 14 )E 2 E 1 E W 1 



Proof. Fix an (m — 2)-length word w consisting only letters 1 and 2. The column 
vector (1, 1, 3) corresponds to the energy measure v 2 obtained from only h 2 . By 
modulo a constant, we have, as shown in the figures of F W SG, the values of h 2 and 
energies on the cells of v 2 . Note that < a < b. 



5a-2b 
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A simple manipulation on matrices shows that 

- I — I .14 )E-,E I E„. [ 1 

3 

15 ( 1 1 1 ) E 2 E 2 E 1 E W \ 1 j - ( 1 1 1 ) E 2 E 1 E W ( 1 



jj (15u 2 (F w F 1 F 2 F 2 SG) - v 2 {F w F x F 2 SG)) . 



Similarly, 

-1 - i 14 ) E\ E 2 E ir \ 1 

3 

15 ( 1 1 1 ) E 2 E X E 2 E W \ 1 j - ( 1 1 1 ) E 1 E 2 E W ^ 1 

= \ {\hv 2 {F w F 2 F x F 2 SG) - v 2 {F w F 2 F x SG)) . 
Whence, from the calculation in the figures, 

-1 -1 14 )E 1 E 2 E W | 1 j - ( -1 -1 14 )E 2 E X E W I 1 



4 /5 
625 U 



(6053a 2 + 852a6 + 1446 2 ) > 0. 



□ 



What we have proven in the above lemmas is that, for two level-m cells F W SG 
and F W /SG whose words consist only of the letters 1 and 2, we have v 2 {F w SG) < 
v 2 {F w i SG), provided that F W SG is to the left of F w >, in our usual orientation. It 
also follows immediately that for those word w of length m, 

15 99 5 /3 s " 

c w > Ci... 1 



5 m 4 • 15 m 2 

which implies (I)' 1 "' c w is bounded away from zero, for every such 
We now turn to an upper bound for f|) w H-E^H, where ||- 1| is the operator 1— norm. 



IV. 



Lemma 5.3. For any harmonic function f and word w with \w\ = m, 

'3 N ™ 



£fM) <2 0sc(f,SG) 2 
Proof. By [g], 

Oscl'/'. F...SG) < ( 



Osc(f,F w SG)< I |) Osc(f,SG). 



18 



RENEE BELL, CHING-WEI HO, AND ROBERT S. STRICHARTZ 



Since for fixed oscillation, the greatest energy is attained by a function which is 
symmetric on F W SG, which has energy (|) m 2 (Osc(/, F W SG) 2 , we have 



£fJ/)<20sc(/,SG) 5 



Lemma 5.4. There exists a constant C, independent of m, such that 

'3" ~ 



□ 



\EJ\<C 



for all \w\ = m 



Proof. Consider || (x, y, z)\\ = 1. If we view x, y, z as the energies of the level 1 cells, 
we know that (a;, y, z) corresponds to an energy measure oo^o + «i^i + «2^2- Since 
(x,y,z) and (00,01,02) differ only by a linear transformation, ||(do, ffli, C2.2) || < B 
for some constant B for all such triples (00,01,02). Let a word w be fixed, with 
\w\ = m. Since from Lemma 5.3 we know that Vj(F wi SG) < 2 (|) m+1 , we have 




J2 J2 a ^(F wi SG) 



<EEI^II^( F ™ 5G ) 

i=o j=a 



<EEI^I' 2 

i=0 j=0 

m+1 



m+1 



< 6B 

Taking supremum over all \\{x,y, z)\\ — 1 and letting C = 6B (|) give 

\E,„\\ < C 



□ 



We now give the main result of this section, the continuity on the three edges of 
a cell. 

Theorem 5.5. Given an energy measure Vf g and a cell F W SG, the restriction of 



the derivative 



dv s, s 

du 



on the three edges restricted to the vertices V* is continuous. 



Proof. Since every point x of V* that is on the edge of a cell is the junction point of 
two cells (boundary points excepted) it suffices to show that the restriction of the 
derivative dv ^ v a on the edge F w [qi,q 2 ] is continuous at F w q 2 , where [01,02] is the 
edge connecting qi and q 2 . Continuity on the other side of x follows from symmetry 
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and arbitrary choice of cell. To prove this, we will prove that 



( -1 -1 14 )E w ,A n 



x 

y 



( - 1 -1 I 4 ) 



lim 




dv f,g 
dv 



{F w q 2 ) 



( -1 -1 U)E w ,A n y/3 j (-1 -1 14)^/3 
uniformly in all finite w' made of 1, 2 where 



Vf.g{F w0 SG) 

v f , g (F wl SG) 
vj, a (F wi SG) 




v{F w0 SG) 
v(F wl SG) 
v{F w2 SG) 



This is enough, since the set {F w >F n }, where w' ranges over all finite words 
consisting only of the letters 1 and 2, is exactly the set of points in the graph V* 
that are in the interval F w F£[qi, q^]- 



Fix a w' with 



and let e > 0. We consider 



( -1 -1 14 ) E W ,A. 




( -1 -1 14 )E w ,A n (3 




( -1 -1 U)E w ,D n [ /3 +(-1 -1 U)E W ,A\ p 

7 / \ 7 



where D n is the "error matrix" of A n . 
By Lemma 5.1 and 5.2, we can see that 



(§) (-1 -1 14 )E W ,A^ p j > ^(14 7 -a-/3)>0 



for all those w'. We also have 

-1 -1 14 ) E W ,A I y 



max 



<||( -1 -1 14 )||||^||P|| max 



< 



5 J C IK -1 -1 14 )|||L4|| ( /3 
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for sonic constant C independent of m by Lemma 5.4. We can choose n large 
enough so that D n satisfies, with lemma 5.4 and that A n — > A in the operator 
1-norm, that 



3 

< 5 



x 



a 



umx<(-l -1 lA)E w ,D n \ -1 U)E w ,D n yf3 



where 5 is chosen small enough that for all 771 , rj 2 < S 



r?i + (!) m (-l -1 U)E W ,A\ y (-1-1 U)E W ,A\ y 



a \ la 

% + (|) m (-l -1 U)E W ,A\ /3 (-1-1 U)E w ,Ai p 

7 / V 7 

which is possible by continuity of the function (x,y) i->- Whence, 



( -1 -1 14 )E w ,A n I y \ ( -1 -1 14 ) [ y 



< e. 



< £ 



a \ I a 

( -1 -1 14 )E w ,A n I ( -1 -1 14 ) i /3 

7 / V 7 

The limit is uniform on all w' consisting of only letters 1 and 2 by the equicontinuity 
of the functions {(x, y) ^ |±f } at (0,0) for < a, & < M, < m < b. □ 



1 

The above idea uses only the lower bound for (|) m ( — 1 —1 14) E w 1 

\ 3 

and upper bound for (I)™ ||-E\o||- 



6. Average Value of the Radon-Nikodym Derivative 

We have shown how to calculate the values of the derivative on the vertices in 
K; moreover, we have continuity on the edges of every triangle. We would like to 
relate the average of derivatives on the whole cell to a weighted average on three 
vertices of the cell; that is, to find, on a cell C, positive real numbers {&j}j =0 wnose 
sum is 1 satisfying 



to) 



for every energy measure of harmonic functions Vj g . In this section, we will show 
that the coefficients bj exist and discover some characterizations of them. 



ENERGY MEASURES OF HARMONIC FUNCTIONS ON THE SIERPINSKI GASKET 21 



6.1. Existence. Recall from Theorem |2.1[ that we have 
f o \ / f o 



M \ v x o F^ 1 + Mi i/i o Ff 1 + M 



^ 2 / \ Vl 





and 

v 2 (FjA) 

Letting M w = M Wl M W2 ■ ■ ■ M Wm we see that 

MF W A) \ 
vi{F w A) =M U 
v 2 (F w A) J 

Theorem 6.1. On each cell C — F W SG, there exist unique (bo,bi,b 2 ) & K 3 with 
bo + b\ + b 2 = 1 such that 

2 




Avg f 



dv f,g 



* ' " ^ j=0 

for all energy measure (not necessary positive) Vf g on SG. Explicitly, 

1 1 sum of column j of M w 



bj - — 

6 2 sum of all entries of M u 



Proof. We first introduce some notation. We let 




, Z = Zq + Z\ + z 2 . 



Now suppose Vfg = a^v^ + ct\V\ + a 2 v 2 . Then we have 

Vf, g {F w SG) = 2 ( a ai a 2 ) z, v{F w SG) = 2 ( 1 1 l)z 

so the average value of the Radon-Nikodym derivative with respect to vj g on the 
cell F W SG is 

If w sg ^T dv = v f ,g{F w SG) = ( a oi a 2 ) z 
If w sg dv v{F w SG) (1 1 1 ) z 

We also note that 

dvj^ v Lg (F w F™SG) 
dv { wQi) ~ m lm io v {F w FySG) 

( pfa) 

( a «i a 2 )M„ | p^(^) 
= lim 



dv Vi3) 



111) M, 




and we know that rf y- g (F w q ) is equal to | if j = fc and \ otherwise. 
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So 

( a ui «2 ) M w I ^ 



dv 



(F w q ) = 



dv 



(1 1 1 ) M u 

( a ai a 2 ) (z + \z) 
(1 1 1 ) (*o + \z) 
( I 

( a ct\ a 2 ) M w I ^ 



(F wqi ) = 



dv 



(1 1 1 ) M u 

( a ai a 2 ) (zi + |z) 

( 1 1 1 ) (*i + P) 
/ l 

( Qf ai a 2 ) M, 



(i^<72) = 



(i i i ) m w y I 

( a ai Q!2 ) (z 2 + \z) 

(i i P) 

Now we suppose that we can write the average of the derivative on the whole cell 
as a weighted average of the value of the derivative on the boundary points; that 
is, we write 



du f>g _ ( a ai a 2 )z 

AV 6F„SG" 



dv (ill 



Then it is implied that 

2 



Z 



b d ^{F wqo ) + b^{F wqi ) + b 2 d ^{F wq2 ) 
i \ ( {zj + \z) 



\ - (jj +_|f) = z 

U( 1 1 1 )fe + P) ( 1 1 1 )*' 



We also see that multiplying both sides of this equation by ( 1 1 1 ) gives 
b Q + b 1 +b 2 = 1. 
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Now we note that the vectors {zj + \z}j are linearly independent since 

M w ^ | j ,M W ^ | j ,M U 

are so by the fact that the matrix M w is invertible, so we have the existence and 
uniqueness of the bj . We now calculate some simple expressions for the bj . 

Now we let 

d 



3 (1 1 

Then we can rewrite the previous equation as 
This means that 



1 



(cfoMo + rfi&i^i + d 2 b 2 z 2 ) 



-j— — / - \ (dobo + dxh + d 2 b 2 ) 

{ 1 1 1 j z o 



so 



d di d 2 



where D is determined by 



1 = b + h + b 2 

bj=D^( 111) {^z+\zi 

1 = D(( 111)^) 
_ ( 1 1 1 )(iz+lz J ) 
(111)- ■ 

Noting that ( 1 1 1 ) Zj is the sum of the entries in column j of M m and 
( 1 1 1 ) z is the sum of all entries of M w . Then 



1 1 sum of column j of M w 

bj = — I 

6 2 sum of all entries of M u 



□ 



Later we will show that b'jS lie between and |; thus, the above theorem shows 
that we can relate the average value of the derivative on the whole cell to the 
weighted average value of the derivative at the vertices of the cell(one set of weights 
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works for all energy measures). Since the 6' s depend on the word w, from now on 



we will denote b j on F W SG as bj 



(w) 



6.2. Characterization. We now have coefficients for the cell F W SG to cal- 
culate the average of derivatives on the whole cell with the values on vertices of 
the cell. Next, we will give two ways to calculate the 6's. One way is to calculate 
recursively from cells to subcells; the other way is to calculate through the Kusuoka 
measure. 

The following theorem gives the first way to calculate all the 6's by calculating 
them recursively. 

Theorem 6.2. Using the previous notations, we have the following relationship 
from cell to subcells: 



96 



13i£ x., 
r »- (w > +2b[ w > 



(t " ) + 13bS™ ) +4 1 



,(w2) _ 26^-b^+2j| 
On - -(»),,(»)_ - ~.<. 



b\ 



(wO) 



b[ W2) = 



2b^ w> +2b\ u 




13b ( m) +b[ u 


>+6<»' 


96^™ 




&£" ) +13&j" 




-b^+2b^ 





'+6^+136^' 

Proof. We know that, on F W SG, 



l(»0) 



6^ 2) = 





+26!™' 


13fc<,™'+6< 1 
-6^ ) +2bf' 




'+24™' 


b<,™' + 136< 1 


"+6<<"' 


9b<™ 


) 



6<™'+6< ra ' + 13&! i ' 



(™) 



1 1 

6 + 2, 



>0 



(to) 



+ c + c. 



(«0 

2 



where 



» 



(«0 _|_ c («0 



+ 4 ,u) ) = 4 u ' ) 



= ( 1 1 1 ) M w . Put = 2(6 
Since 



J ) . We see that 



-o c i 



4"° ) M = ^ ( 94 w) + 2c^ ) + 2c> ) 2c> ) - c! 



(«,) 

2 



2c ( ™ ] ~ c^ ] 



we have 



and hence 



2c 



1 1 94 w) + 2c^ ) 
6 + 2 9c (-) + 3c («0 + 3{ >) 



°0 



Similarly, we have 



and 



b\ 



(wO) 



1 1 9a + 2ai + 2a 2 
6 + 2 9a + 3ai + 3a 2 ' 



1 1 2ai — a 2 
6 + 2 9a + 3ai + 3a 2 



M) = 1 1_ 
2 6 2 9a 



2a 2 



Oi 



3ai + 3a 2 



With some computation from the definition of a's and using the fact that J2j = 
1, we can relate 6^ and b[ w0 ^ explicitly: 



°0 



96^ } 


I3b { ™ } + 




_ 2b ( w) + 2b[ w) 




m Q w) + b[ w) 
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,(w0) _ Zo Q — Q\ + M 2 



+ b[ w) + b { ™ ] ' 



I •->/ <»•' i ;,' ' 

Observe that multiplying different Mj 's only differs by some permutation of sub- 
script and entries, we have the rest of results. □ 



The next theorem presents the second way to calculate b^ w \ 

Theorem 6.3. The distance ofbj from 1/3 is proportional to how "skewed" the 
Kusuoka measure is on the cell F w FjSG relative to F W SG. That is, 

■ [b{w) _l\ = l{v(F w F j SG)_l 



4 V y{F w SG) 



Proof. We first simplify the term 

( 1 1 1 ) M w I I (1 





( 1 1 1 ) M w 1 ( 1 



Note that since 2 = v (SG) = V\{SG) = ^(SG), the denominator is just 
v(F w SG). Next, we see that 



(1 1 l)M w \ 




4/5 

1 11) _U„ i 





.(!)((. i dm. (| 

h - (v(F w F SG) - l -v{F w SG) 



2/5 

(1 1 l)M w \ 2/5 

2/5 



where the last equality comes from noting that 6/5 = v (F SG) and 2/5 — vi(F SG) 
v%{F SG). So we have that 



,( W ) v(F w F 3 SG)~lv(F w SG) \ = 5 f -I v(F w FjSG) \ 

3 4 v 15 v{F w SG) ) 4 V 15 v{F w SG) ) 

which gives 

1 (,(w) _ 1\ = 1 ( u(F w FjSG) _ 1\ 
5 V j 3/ 4 v v(F w SG) 3 J ' 

□ 
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As a consequence of the theorem, the bounds for (b^ — |) also give bounds for 
how skewed the Kusuoka measure is in a cell. We now provide some bounds for the 



coefficients b 



(w) 



Theorem 6.4. a) The infimum of {b^} over all words w is 0. 
b) The supremum of {b^} over all words w is |. 

Proof. We first prove a). We know that 
1 



15 



14 



-1 )E U 



1 

15 



lim 



v{F w F™) > 



and the inequality is sharp over all w because every energy measure is continuous. 
So 



-(14 -1 -1 )E U 



= ( 1 1 1 ) (Eo ~ ^I)E U 



= ( 1 1 1 ) E E U 
> 



15 ( 1 1 1 ) E « 



and hence 
(111) E E U 

which implies 



(1 1 l)E u 



>- 15 ( 1 1 l)E u 



v(F w FySG) - \v(F w SG) > _2_ 
v{F w SG) ~ 15 



and so finally we get that 



5^2 v {F w F SG) - \v(F w SG) 
u{F w SG) 



4 V 15 



> 0. 



Since the inequality is sharp, zero is the infimum, completing a). 
To prove b), we see that 



( -2 3 3 ) E u 



1 



( -1 -1 14 )E 1 E W 2 = - lim 



HF w ,iF™) > 



and this inequality is sharp. So 
5iy{E w F SG)-3v{F w SG) = ( 1 1 1 ) (5E -3I)E U 
Whence 



(2-3-3)^1 2 )<0. 



b -{v(F w F SG) - \v{F w SG)) < v{F w SG) 
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which implies 

5 v{F w F SG) - \v{F w SG) 
2 v{F w SG) 

and so finally we have 

{w) = 1 5 is(F w F SG) - \u{F w SG) 2 
6 4 v{F w SG) "3' 

Since the inequality is sharp, this is the supremum. We have b). □ 

The previous results concern the bounds for a single value .It is natural to 

consider (b^"\b^"\b^) as a vector in M 3 , lying in the plane {(x,y, z) : x + y + z = 
1}. Surprisingly, these vectors all lie in the disk centered at (|, |, |) of radius 

Theorem 6.5. For all words w, 

2 

< 



6 



and i/iis inequality is sharp. 

Proof. Let ( 1 1 1 ) M ffl = ( 4"° ) • Then we can have 



and so 



4 

So we need to show that 




(to) 



cq T4 -r '--2 



1e 



4 ~ 4 

for which it suffices to show that 

„{w) (w) (w) (w) (w) (w) n 

which we will prove by induction on the length of w. When \w\ = 1, c^c^ + 

_»„(«>) (w) (w) _ 13 i 13 i i i „ 
c l c 2 + c c 2 — 15 15 + 15 15 + 15 15 ^ U " 

Now suppose it is true for w with \w\ = m. We see that 

( ci w0) c^ 0) 4 w0) ) = ( <>> <>> ) M 

= ^ ( 9c { w) + 2c^ ] + 2c { 2 w) 2c{ w) - c { ™ ] 2c 2 w) - c[ w) ) 
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and so 



15 2 (4 w0) C ^ 0) + c^ 0) c^ 0) + c^ 0) 4 w0) ) = 9c ( ^(c ( ^ + 4 W) ) + 2{c^ ) f + 2(4" } ) 2 

+ Ac^c^ + hc^c^ - 2(c{ w) ) 2 - 2(c { 2 w) ) 2 

= 9(c^ ) c^ ) + c{ w) c 2 w) + c { w) 4 w) ) 
> 0. 

As for the sharpness of the inequality, consider a word w of length m consisting 
of only the letter 0. We have that 



°0 — O Olo2m. I 1\>°1 - °2 - 



1 



3(3 2m + 1)' 1 2 3(3 2m + 1) 6 
and so as m oo we have lim m — \) 2 = \. □ 



Remark 6.6. In Theorem 6.2 we see that from F W SG to F w qSG, the properties of 



bj were translated to the properties of the rational transformation Aq 

9ao+2ai+2a2 
9ao+3ai+3a2 



/ «0 




Ol 


l-l 


V «2 , 





9ao+3ai+3a2 
2a2 — ai 



9a +3ai+3a 2 

and we have Ai and for the maps from F W SG to F w iSG and F W 2SG. Thus it 
is worthwhile to look at the properties of aj. First, it is easy to see that Yl a j = 1- 
Also, the sum of squares of aj is strictly less than 1. To see this, using aj — 

2 fb^ - i), we have 



1\ 2 1 

<6 



Since ^2a,j — 1, expanding the above gives 



2 

3/ <3- 



Finally, in terms of the ratios of Kusuoka measure, we have 

u(F wj SG) _ 2 / _ 1 



v{F w SG) 5 V 2 



7. Energy Distribution 



Let &("> = (6^,6^,6^). We have seen in Theorem 6.3 how the vector b^ 



relates to the splitting of Kusuoka measure in the cell F W SG when it subdivides 



into three subcells at the next level. Theorem 6.2 gives a recursive algorithm to 
compute the b^ w > vectors. In this section we investigate the distribution of the col- 
lection of all vectors as w varies over all words of fixed length m. We may write 
b w0 = B {b w ), M™ 1 ' = Bx(bW) and b^ = B 2 (b^) where the Bj are rational 
map given in Theorem |6.2| We may regard {B , B±, B 2 } as an iterated function 
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system acting on the open disk described in Theorem |6. 5 1 



(6.0) r sin(4>) (2.0/3 0)+sqrt(2.0) sqrtp.O) r sin(<B (2.0/3.0J+1 .0;3.0)/(sqrt(2.0) r cos(4>) 6.0+sqrt(B.O) r sin( 




Figure 7.2. The map B 



In Figure 7.2 we illustrate the geometric structure of Bq by showing the image 
under Bq of a uniform polar coordinate grid on the disk. B\ and B2 are similar, 
only rotated through angles ^ and = ^ L . 



To obtain {b {w *> : \w\ = 1 
vector (~, ~, |). In Figure 
this distribution of points in t 



} we apply all 3 m iterates of {B , Bi,B 2 } to the initial 
we show the result for m = 12. To better understand 
re disk we examine its angular and radial projections. 



7.3 



Introduce polar coordinates (r, 9) for the disk with r < (the origin corresponds 
to the vector b = (|, |, |)). The angular districution at level m is 

P m (A) = 3" m #{points with 9 e A}. 
By symmetry it suffices to understand the distribution for A C [0, ^]. In Figure 



we show histograms (100 slices) of P m on [0, =jf\ for m — 11 and m = 13. 



Conjecture 7.1. Let P m be defined as above. The limit of P m as m — > 00 is an 
absolutely continuous (Bq, B\, B2) -invariant measure. 

Similarly, the radial distribution at level m is 

Q m (A) = 3" m #{points with r e A}. 
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Figure 7.3. Points of {&("') : \w\ = 14} 



In Figure 
m = 14. 



7.5 



we show histograms (300 bins) of Q m on [0, ^] 



for m — 10 and 



Conjecture 7.2. Let Q m be defined as above. The limit of Q m as m — > oo is the 
delta measure on the boundary r = . 

We can say more about the mappings. Recall that we have defined 



B (x,y,z) 



9x 2x + 2y - z 2x - y + 2z 

13x + y + z' l'Sx + y + z' 13x + y + z 



by Theorem 6.2 for each word u>, we have B (bW) = b^ w0 l 

We first introduce the "polar coordinates" on the disk centered at (|, |, |) on the 
plane {x + y + z = 1}. We know that (cos(6>+ |), sm(0 + ^), 0) is a parametrization 
of the unit circle lying on the x — y plane. Put 



.4 



1 V2 ^6 

o - 2 



i 

t 



\/6 V3 



which maps the unit circle lying on the x — y plane to the unit circle lying on the 
plane {x + y + z — 0}. It follows that the parametrization of the circle of radius r 
centered at (I, |, I) on the plane {x + y + z = 1} is 



(7.1) 



.4 



rcos(6> + |) 
rsin(6>+ f) 




v'TT 



cos ( 



: sin u — 



■- cos I 



^isin0-^cos( 
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(a) Level 11 (b) Level 13 

FIGURE 7.4. Angular Distribution 



50 100 150 200 250 300 



50 100 150 200 250 300 350 



(a) Level 10 (b) Level 14 

Figure 7.5. Radial Distribution 



With this parametrization, the image of radius r under Bq is 

cos 9-\-S \ 



(7.2) 



cos 04-8 
4\/6rcos0+5 

cos v— 4 
4^1-008 0+5 
-9V2r sin 6-5V6r cos 0-4 
4\/6rcos0+5 



Thus if we take r = , we see that the image of the boundary circle is 
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Each Bj maps the disk onto itself injectively. We first show that the mappings 



map the disk into itself. By (7.1) and (7.2), we need to see, by considering B , 
whether to each pair (r, 9), there exists a unique solution (7, a) satisfying 

/ 10 V6r cos \ / r> /c 

1 -^—^ \ I 27V6cosa 

37V2 sin a — 7v6 cos 9 



V 



4vl3rcos6l+5 

cos 6 — 4 
4%/fjr cos 9+5 
-9\/2r sin d-bV&r cos 8-4 
4V6rcos0+5 



\ — 37V / 2sina — 7-^/6 cos 



The above reduces to 
(7.4) 

In particular, 
(7.5) 



7 cos a 



7 sin a 



1 5\/6r 



V6 4\/6r cos 6> + 5 
3r sin 

4v6r cos + 5 



5) 2 



so if r < 4g, then 7 < and there is a solution (7,0) satisfying (7.4). But by 
(7.4), we see that actually (r, 6) is also completely determined by (7, a); it follows 



that the Bj maps the disk onto itself injectively. 



When we restrict the map on the boundary circle, we see from (7.5) that the 
circle is mapped onto itself and using the notations as before, 

5 cos 9 + 4 
4 cos 9 + 5 
3 sin 9 



cos a 



(7.6) 



sin a 



4 cos 9 + 5 

These show that a corresponds to exactly one 9; the map gj, referring to the restric- 
tion of Bj to the boundary circle, given by gj{9) = a is well-defined. Differentiate 

It follows that 



dg a (e) 



both equations in (7.6), we have — ^g- 

3dt 



9o(0) = 



4 cos t 



4 cos 0+5 ■ 



= 2 tan 



1 



■ tan ■ 



and 



% V) 



2 tan 



1 (3tan- 
V 2 



Since g\ and gi differ from go a rotation, we have 
g x (e) = 2tan~ 1 

gi 1 (a) = 2tan _1 



1 

- tan 
3 



3 tan 



g 2 (9) = 2 tan 1 ( - tan 



g 2 1 (a) = 2 tan 1 ^3 tan 
We have the following theorem. 



2tt 

T 

2?r 

T 

2tt 

T 

2tt 

T' 
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Theorem 7.3. The mappings Bj (j = 0, 1, 2) map the disk U centered at (|, g, |) 
of radius ^= onto iiseZ/ injectively. Further, 

(i) £/ie mappings map the boundary circle C centered at (™, |, ~) of radius 
onto itself; and 



(ii) If is the function in (7.11 which associates each (^g>$) £ R to C, t/iere 



are differ entiable maps gj satisfying Bj(^{^,8)) = ^(^Tg, 9j(@)) given by 

( 1 
ffo (0) = 2tan~ 1 ( - tan- 

^) = 2tan- QtanQ-|))+f 

ff2 W = 2tan- Qt-(^ + |))-f 

(iii) -Eac/i -B^ /las exactly two fixed points on the circle. The fixed points are, in 
"polar coordinates", *(^,0), *(^,7r) for B , f ), ^) /or 

Proof. We have already proved (i) and (ii). For (iii), we only need to prove for So 
and the other two follow. By (ii), 9(A^,0), \f(-^,7r) are the only fixed points on 
the boundary circle. So it suffices to show that there is no fixed point in the open 
disk. Suppose, on the contrary, we have such a fixed point, say ^(r, 9). Then we 
must have _ 

1 lO-v/frr cos 9 + 8 2r 

= = — = cos 9 

6 4yV cos 6 + 5 

which shows, since r < -^=, 

cos 2 9 = -L > i 
6r^ 

but this contradicts the fact that | cos£?| < 1. □ 
By (ii), We can regard Bj and gj on the boundary circle as the same map under 



the parametrization (7.1). 



If we directly view the maps Bj as maps on the unit disk, we have an easy 
description on how they act on the unit disk. 

Theorem 7.4. Using the notation (r,0) and (7, a) as before. Let (x,y) = 
(\/6r cos#, y/6r sin#) and (x',y r ) = (v^T cos 1/67 sin #) be on the unit disk. Then 
Bo maps every vertical line to a vertical line and every point (x,y) with x > — \ 
closer the boundary circle. In other words, if we denote T to be the triangle with 

vertices (1,0), (—5, 2 )> ( — 5' ~2 )> ^ ien ever V point in T is mapped closer to the 
boundary by all of the three mappings Bj and every point outside T is mapped closer 
to the boundary by two of the mappings Bj . 

Proof. The map Bq on the unit disk is of the form 

, 5x + 4 3y . 



B {x,y) = Or/, t/) = ( 



4x + 5 ' 4ir + 5 ' 
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Fix a constant c G [—1,1]. If x = c, then x' = ffzf ; so i?g maps every vertical line 
to a vertical line. Also, by (7.5), we have 



2 1 _ 9 / 2 1 



7 "= r^a l r - 



6 (4x + 5) 5 



thus if x > — = , ,„ ? c . 2 < 1 and hence 

2 ' (4a;-|-5) 



2 1 9 ( 2 1 

7 < r 

' 6 (4x + 5) 2 V 6 



The last assertion comes from the symmetry of Bj, since all of these mappings 

,,,,,, - r\r\\\r 

3 " " 3 



differ by a rotation of and — only. □ 



By Theorem |7.3[ we can say more about the measure mentioned in Conjecture 



7.1 We call the measure dX(t) = f(t)dt for some /, X[a, b] — X(Bj[a,b]) for all 



j = 0, 1, 2. Recall that it is (B , Bi, B 2 ) -invariant, hence 



[ /= / f(t)dt 

J a J q- \a,b] 



l [a,b] 

It follows that the function / can be characterized by one satisfying the relation 



/(*) = 

j = Q »J 



If we consider gj separately, we have 



1 4 /„._!/„. t\\ 5 



= - cos 2 tan 3 tan 



5o(5o _1 W) 3 V V 2JJ 3 5-4cost' 

1 _ 3 
flifor'W) ~ 5-4cos(i-f) ; 
1 3 



fliGfeH*)) 5-4cos(i+f)" 

Finally, we give further experimental evidence for the conjectures. Conjecture 

7.2 states that the radial distribution is the delta measure at r = the angular 
distribution on the whole disk should be the angular distribution on the boundary 
circle only. Pick three points (0,1), (—1,0), , — ^75) from the three arcs men- 
tioned in Theorem 7.4 respectively and iterates them 14 times by all of the three 
maps. If we cut the whole circle into 800 slices and plot the histogram on one sixth 
of circle the number of points fallen into the corresponding bins as in Figure [TU] we 



see that the picture is extremely close to the histogram of level 13 in Figure 7.4b| 
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Figure 7.6. Angular Distribution on Boundary 
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